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Abstract 

Background: Di-(2-ethylhexyl)-phthalate (DEHP) is a commonly used plasticizer in polyvinylchloride (PVC) 
formulations and a potentially non-genotoxic carcinogen. The aim of this study was to identify genes whose level 
of expression is altered by DEHP by using a global wide-genome approach in Syrian hamster embryo (SHE) cells, a 
model similar to human cells regarding their responses to this type of carcinogen. With mRNA Differential Display 
(DD), we analysed the transcriptional regulation of SHE cells exposed to 0, 12.5, 25 and 50 liM of DEHP for 24 hrs, 
conditions which induced neoplastic transformation of these cells. A real-time quantitative polymerase chain 
reaction (qPCR) was used to confirm differential expression of genes identified by DD. 

Results: Gene expression profiling showed 178 differentially-expressed fragments corresponding to 122 genes after 
tblastx comparisons, 79 up-regulated and 43 down-regulated. The genes of interest were involved in many 
biological pathways, including signal transduction, regulation of the cytoskeleton, xenobiotic metabolism, 
apoptosis, lipidogenesis, protein conformation, transport and cell cycle. We then focused particularly on genes 
involved in the regulation of the cytoskeleton, one of the processes occurring during carcinogenesis and in the 
early steps of neoplastic transformation. Twenty one cytoskeleton-related genes were studied by qPCR. The down- 
regulated genes were involved in focal adhesion or cell junction. The up-regulated genes were involved in the 
regulation of the actin cytoskeleton and this would suggest a role of cellular plasticity in the mechanism of 
chemical carcinogenesis. The gene expression changes identified in the present study were PPAR-independent. 

Conclusion: This study identified a set of genes whose expression is altered by DEHP exposure in mammalian 
embryo cells. This is the first study that elucidates the genomic changes of DEHP involved in the organization of 
the cytoskeleton. The latter genes may be candidates as biomarkers predictive of early events in the multistep 
carcinogenic process. 



Background 

Di-(2-ethylhexyl)-phthalate (DEHP) is a commonly used 
plasticizer in polyvinylchloride (PVC) formulations 
which have a number of applications, especially in food 
packaging, medical devices or cosmetics. Phthalates are 
not chemically bound to PVC and can migrate from 
PVC-containing products to the environment, resulting 
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in significant environmental contamination and human 
exposure [1,2]. 

DEHP experiments have revealed toxicities including 
carcinogenesis and endocrine-disrupting effects, but no 
genotoxicity has been recorded [3]. DEHP is capable of 
disturbing the reproductive process by mimicking or 
antagonizing steroid hormone action [4] and its effects 
on testosterone, luteinizing hormone or estrogen-like 
activity have been reported [5,6]. DEHP has been shown 
to decrease free testosterone levels in humans after 
occupational exposure [7] and thyroid hormone levels in 
adult men otherwise exposed [8]. 
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DEHP has been classified as a peroxisomal proliferator 
and as a non-genotoxic carcinogen in animals [9]. 
Experimental studies using rodents and in vitro assays 
showed that DEHP and its active metabolite MEHP 
(mono-(2-ethylhexyl)-phthalate) can interact with 
nuclear receptors like PPARa [10] or PPARy [11]. Oxi- 
dative stress, as a result of peroxisome proliferation, and 
DNA damage have been described in the human pros- 
tate adenocarcinoma cell line LNCaP [12,13] and the 
mouse Leydig tumor cell line MA-10 [14] exposed to 
high concentrations of DEHP (3 mM). Peroxisome pro- 
liferation is one of the mechanisms that produce liver 
tumors in rats or mice, but this mechanism was not 
judged to be relevant in humans [15]. The liver is not 
the sole target for DEHP carcinogenicity: testicular 
tumors [16] and pancreatic acinar adenomas have also 
been reported [17]. Other studies have pointed out that 
peroxisome proliferation is not a necessarily pathway in 
the carcinogenicity of DEHP [18] and more liver tumors 
occurred in PPARa-null mice than in wild type animals 
[19]. Transcriptional changes independent of PPARa 
were also found in rats and mice exposed to DEHP [20]. 
Several non-PPARa mechanisms were addressed: activa- 
tion of p38 mitogen-activated protein kinase not 
involved in peroxisome proliferations [21]; stimulation 
of growth regulatory pathways, mitogen-activated pro- 
tein kinase, extracellular signal- regulated kinase and p38 
phosphorylation [22]. Other mechanisms related to non- 
genotoxic carcinogenicity, like inhibition of gap junc- 
tional intercellular communication [23] or inhibition of 
apoptosis, were reported. Apoptosis was shown to be 
suppressed by DEHP through different pathways. An 
interference with the cytokine TGF-pi (transforming 
growth factor-pi) [24] or with TNF-a (tumor necrosis 
factor-a) has been described [25]. An increased level of 
Bcl-2 and negative regulation of c-Myc expression has 
been related to inhibition of apoptosis in Syrian hamster 
embryo cells treated with 50 [iM of DEHP [26]. 

Several authors have demonstrated that DEHP and its 
active metabolite MEHP induce morphological transfor- 
mation of SHE cells [27-29], indicating the carcinogenic 
potency of the two chemicals. Although phthalate toxi- 
city has been extensively investigated over the past 10 
years, the mechanisms of DEHP carcinogenicity have 
not been elucidated. It was recently stated by the Inter- 
national Agency for Research on Cancer (IARC) that 
PPAR-independent mechanisms of DEHP carcinogenesis 
are necessary to be studied [http://monographs.iarc.fr/ 
ENG/Publications/techrep42/TR42-18.pdf]. 

The choice of cellular models and methodologies is 
critical to the study of the phenomenon of carcinogen- 
esis. Syrian hamster embryo cells are a relevant model 
for mechanistic studies of chemical carcinogenicity. SHE 
cells, unlike mouse and rat cells, are less responsive to 



peroxisomal proliferation and, in this respect, more 
similar to human cells. SHE cells are normal, diploid, 
genetically stable and primary cells which are metaboli- 
cally competent for procarcinogen activation. Therefore 
they are used to study mechanisms of in vitro carcino- 
genesis [30]. SHE cells are obtained from embryos after 
removal of the differentiated tissues, and the population 
is mainly composed of epithelial and fibroblastic cells 
[31]. SHE cells from colonies having been morphologi- 
cally transformed after short exposure to chemical carci- 
nogens induced tumours when transplanted back into 
hamsters [32]. This validated the model and the cell 
transformation criteria for in vitro carcinogenicity. 
Recently, the SHE cell transformation assay has been 
recommended by OECD in 2007 as in vitro method of 
screening chemical carcinogens on the basis of its per- 
formances to detect non-genotoxic as well as genotoxic 
carcinogens [http://www.oecd.Org/dataoecd/56/5/ 
37863750.pdf]. 

The aim of this work was to use a global transcrip- 
tomic approach to understand the molecular mechan- 
isms of cell transformation induced by DEHP in SHE 
cells. The objectives were to identify changes in gene 
expression occurring in the early steps of cell transfor- 
mation as well as pathway disturbances that may trigger 
a carcinogenic process. A characterization of the genes 
expressed in SHE cells at DEHP concentrations inducing 
cell transformation may give information on PPAR-inde- 
pendent mechanisms and alternative pathways of DEHP 
carcinogenicity. The transcriptomic changes induced by 
DEHP in SHE cells were analyzed in the first hours of 
exposure. We focused secondly on changes of cytoskele- 
ton-related genes underlying morphological transforma- 
tion in SHE cells. Indeed, cell transformation is 
expressed by the alteration of cell morphology, a disor- 
ganized pattern of colony growth and the acquisition of 
anchorage-independent growth which is predictive of 
their ability to induce tumors when injected into syn- 
genic animals [33]. Despite the central role of the actin 
cytoskeleton throughout the life cycle, little is known 
about the gene expression changes involved in deregula- 
tion of its dynamic in the first stages of tumorigenesis. 
Cytoskeleton defects in relation to cancer have been 
mostly studied in the late stages of cell invasion and 
metastasis. 

Differential Display was chosen to identify differen- 
tially-expressed genes in SHE cells and to explore the 
entire genome. The mRNA differential display described 
by Liang and Pardee [34] is a powerful approach for 
transcriptomic analysis. This methodology has become 
popular as a tool for non-model organisms because of 
lack of requirement of previous genomic information 
about the species of interest. As the genome of hamster 
is partly characterized so far, Differential Display 
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appeared quite appropriate to study DEHP dose-depen- 
dent effects in SHE cells. We applied the current metho- 
dology that uses a combination of 3 "anchored" oligo-dT 
primers (to divide the cDNA population in 3 subsets) 
and 80 "arbitrary" primers of 13-mers. The 240 primer 
combination allowed us to obtain a level of 95% gene 
coverage [35]. 

DD was applied to cells exposed for 24 hrs to DEHP. 
Genes corresponding to differentially-expressed frag- 
ments were characterized. Differential Display allowed 
us to screen a set of differentially-expressed fragments 
in treated cells, among which 122 genes were identified 
as targeted by a 24 hr-DEHP exposure. These genes 
were involved in such functions as transcription signal- 
ling pathways, cytoskeleton regulation, apoptosis, meta- 
bolism. As Differential Display is a semi-quantitative 
method, the expression changes of the genes we were 
interested in, were checked by qPCR using hamster spe- 
cific primers. qPCR was applied to RNAs not only from 
24 hr-treated cells, but also from cells treated for 5 hrs 
in order to study the cell response in the meantime. We 
particularly focused on changes of cytoskeleton-related 
genes underlying morphological transformation in SHE 
cells. The objective was to explain from a mechanistic 
point of view the gene expression changes after DEHP 
exposure. To the best of our knowledge, this exercise 
has never been done previously. 

Results 

Identification of DEHP-responsive genes using Differential 
Display 

The Differential Display technique was used to identify 
genes differentially expressed in SHE cells, after 24 hrs 
of treatment with DEHP. An illustration of differen- 
tially-expressed fragments is given in Figure 1 which 
shows gels obtained after the DD protocol and high- 
lights fragments regulated more than 2-fold by DEHP. 
Using 3 anchored primers and 80 arbitrary primers, 178 
differentially expressed fragments were identified (115 
up-regulated and 63 down-regulated). Among these 
transcripts, 141 (79%) showed homology to known 
genes in the RefSeq database (mouse or human) of Gen- 
bank, while 37 (21%) had no homology or homology to 
hypothetical proteins. The sequences of the fragments 
obtained by DD have been deposited in the Genbank 
dbEST database (LIBEST_027390 Syrian Hamster 
Embryo cells library). These 141 fragments corre- 
sponded to 122 genes that are listed in table 1 with 
their accession numbers and the tblastx expected (the 
threshold for tblastx comparison was p < 0.001). These 
genes were classed according to 8 biological functions 
with reference to the GO process database. These func- 
tions included signal transduction and transcription, 
cytoskeleton regulation, xenobiotic metabolism, 
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Figure 1 Representative Differential Display. This figure shows 
results obtained with control (C), vehicle control (VC) and DEHP- 
treated (12.5, 25 and 50 uM) SHE cells mRNA, after 24 hrs of 
exposure. Arrows indicate fragments regulated more than 2-fold 
and used for the next step of analysis. The DD fragments were 
generated using the A-anchored primer (H-dT+A) combined with 
the arbitrary primer H-AP37 (lane a) or H-AP52 (lane b) and the G- 
anchored primer (H-dt+G) combined with the arbitrary primer H- 
AP40 (lane c). After characterization by sequencing, we obtained the 
following genes: al :Rlp37a/a2:Rlp27a/a3:Col 1 al/bl :Cttnbip1/b2: 
Enah/c1: CorolC. 

^ J 

apoptosis, lipidogenesis, protein conformation or trans- 
port and cell cycle. The regulation of the cytoskeleton 
was one of the most impacted pathways. Indeed, 21 
genes involved in this function were differentially 
expressed after DEHP exposure. Ten genes were up- 
regulated, and 11 were down-regulated. 

Transcription and signal transduction is another biolo- 
gical process targeted by DEHP treatment. We found 22 
up-regulated genes, among which 3 were up-regulated 
more than 10-fold (rablb, a Ras oncogen family mem- 
ber, Homeobox 1 and Forkhead PI). Heat- shock 
response related genes (hsp90, hsp70 and hsphl) and the 
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Table 1 List of genes identified by Differential Display after DEHP exposure and classified by major biological 
function, according to GO process 



Id 


Official ^vmhol 


Gene Name 


Accession Nr 


Thla^tv pynprt 


Effect of 






1 -Transcription and Signal Transduction 






AUOU 1 


nspyu 


Heat shock protein 90 kDa protein 


NM_0 10478 


yt-zj 


+ 


MUOUZ 


|_| cn 7fl 
nsp/U 


Heat shock protein 70 kDa protein 


NM_005346 


I t-z^- 


+ 


M jZU 1 


C6bpd 


CCAAT/enhancer binding protein delta 


NM_005195 


qf 1 n 

yC- I U 




a a i m 
1 Uz 


NCDpz 


nuclear cap binding protein subunit 2 


NM_007362 


AP 0 1 
ut-Z I 


+ 


Acini 
MD I U I 


Mottle 


methyltransferase like 5 


NM_029280 


7P 1 n 
/ 1- 1 U 


+ 


AC1 f\") 


Pint 


0-linked N-acetylglucosamine transferase 


NMJ81672 


JL ZO 


+ 


MO 1 U 1 

A6102 


Mr1 \1 
INI I IZ 


nuclear receptor subfamily 112 


NM_003889 


1 F 1 
1 1- I JO 

3E-74 


+ 


A6302 


Mapk15 


nilLUytrN dLLIVdLcU piULclll KIMdSc ID 


NIK/I 177Q99 
I N J VI I / / ;?ZZ 


3E-16 


+ 


A6801 


Rablb 


niclllUcI Ul nrtj UNLUycNc Idlllliy 


MM n^HQRI 


1E-22 


+ + 


A7502 


Mdb1 


iiicuiyi ^_po uiiiuniy uuiiidiii piuicin i 


MM 013^04. 
|\|IVI_U 1 DJy'-t 


2E-16 




C0602 


Gper 


O piULclll LUUjJIcU CollUycll IcLc|JIUI I 


mm nnmQR?m 

I NIVI UU I U.70ZU I 


6E-68 


+ 


C1001 


Hmboxl 


nunicuijux i 


MM CMA^fsl 
I \IIVI_UZ^OO/ 


2E-32 


+ + 


CI 002 


Mapk3 


[IllLUytMl dLUVdlcU piULclll Kllldbc D 


mm nm 1 nQRQi 

I NIVI UU I I \JyOy I 


9E-60 




C1301 


Usp3 


UUiyUILIll pcpLIUdSc D 


mm nn^37 


9E-34 


+ 


C1401 


Psmc5 


26S protease regulatory subunit ATPase 5 


NM_002805 


2E-23 




C4601 


Hsphl 


heat shock 105kDa/1 10kDa protein 


NM_006644 


5E-07 


+ 


C4802 


Irf2 


interferon regulatory factor 2 


NM_008391 


8E-19 


+ 


V»/ IU 1 


Lmo4 


LIM domain only 4 


NM_006769 


7E no 

/t-uy 


+ 




jll IdILL 1 


SWI/SNF related actin dependent regulator of chromatin 


NM_003074 


3F 1C\ 

JL ZU 




G1001 


Lna4 


ATP dependant helicase 4 


NM_001273 


1 E-45 




G1504 


Foxpl 


Forkhead box P1 


NM_032682 


3E-50 


+ + 


VJ 1 / U 1 


L-IcUjI I 


cAMP responsive element binding protein 3-like 1 


NM_052854 


PF 1 Q 

oQ- 1 y 


+ 


vj i yuz 


Pbrm 


Polybromo domain 


NM_001081251 


RF AQ 

on ^-y 


+ 


vjooU 1 


RvfrO 
nXTpz 


relaxin/insulin-like family peptide receptor 2 


NM_1 30806 


1 F 


+ 


VJODU 1 


AKapj 


A kinase anchor protein 5 


NM_004857 


1 f 


+ 


G3801 
G7001 


Eif1 


eukaryotic translation initiation factor 1 


NM_005801 


3E-42 
1E-15 


+ 


G4801 


Spry3 


bpiuuiy iiuniuiuy d 


mm nn^RA-C) 

1 NIVI UUjo4U 


0.0002 


+ 


G6001 


Foxa3 


fr^rl/hioaH hnv A^ 
lUllxllcdU UUA r\D 


MM 00^7^0 
1 NIVI uuozuu 


9E-17 




G6201 


Mett5d1 


rnpthx/ltran^fpr^^p S rnnt^ininn 1 
1 1 icli ly i li cm ioi ci doc -J l.ui i Lai i 1 1 1 ly i 

Z - r\t;y UldUUll Ul LyiUaKcltrlUll 


NM 0/9790 


0.0038 


+ 


A0201 


Cttnbp2 


L-UILdLLIN UlllUlliy piULclll Z 


kin/1 n^n?AQ 

|\|IVI_UjUZ^-y 


0.0001 


+ 


A1901 


Snx6 


suiLiiiy iicxiri o 


KIN/1 091 1AQ 
INIVI_UZ I 


8E-10 




A2401 


Lrrc8a 


IcULIIIc IILII IcpcdL LUI ILdll III iy Or\ 


MM 17777S 

InIVI I / / / ZJ 


5E-32 




A3501 
A3502 
C3501 
G3502 
G3503 


Actb 


beta-actin 


NM_007393 


3E-51 

1E-106 

2E-21 

7E-114 

3E-48 


+ 


A3702 


Collal 


collagen, type I, alpha 1 


NM_000088 


1E-21 


+ 


A3801 


Nrp2 


neuropilin 2 


NM_00 10774 


1E-22 




A5201 


Ctnnbipl 


catenin beta interacting protein 1 


NM_020248 


0.0008 




A5202 


Enah 


enabled homolog 


NM_008680 


2E-09 




A5402 


Kif23 


kinesin family member 23 


NM_024245.4 


1E-25 


+ + 


A7501 


Cdh3 


cadherin 3 


NM_007665 


0.0007 




C4502 


Nid2 


nidogen 2 


NM_008695 


0.006 




C6603 
C6604 
C6605 


Cripl 


cysteine-rich protein 1 


NM_007763 


1E-31 
9E-08 
3E-61 


+ 


G0601 


Thy1 


Thy-1 cell surface antigen 


NM_009382.3 


0.001 
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Table 1 List of genes identified by Differential Display after DEHP exposure and classified by major biological func- 
tion, according to GO process (Continued) 



G0801 


Calml3 


calmodulin-like 3 


NM_027416.3 


0.0004 


+ 


G1102 


Flrt2 


fibronectin leucine rich transmembrane protein 2 


NM_201518 


2E-08 




G1301 


Has2 


hyaluronan synthase 2 


NM_008216 


8E-18 




G1901 


Plekha5 


Pleckstrin homology domain A5 


NM_019012 


3E-06 


+ 


G2301 


Thbsl 


thrombospondin 1 


NM_011580 


3E-20 


- 


G4001 


CorolC 


coronin, actin binding protein 


NM_014325 


4E-07 


+ 


G4802 


Tubb2b 


Tubulin beta 


NM_023716 


0.0001 


+ 


G6901 


Ddkl 


doublecortin-like kinase 1 


NM_019978 


1E-06 


+ 






3-Xenobiotic metabolism 








A0901 


Cyp2e1 


Cytochrome P450 2e1 


NM_021282 


3E-44 


+ 


A2901 


Ephxl 


Epoxide Hydrolase 1 


NM_000120 


4E-21 


+ 


A5701 


Gstpl 


glutathione S-transferase, pi 1 


NM_000852 


1E-24 


+ 


A6601 


Gsttl 


glutathione S-transferase, theta 1 


NM_008185 


1E-44 


_ 


C0101 


Txnrdl 


Thioredoxin reductase 1 


NM_001 042523 


1E-20 


+ 


C1501 








1E-14 




G1502 








3E-79 




C2901 


Cypl b1 


Cytochrome P450 1b1 


NM_009994 


0.0 


+ 


C3101 


Gstm5 


glutathione S-transferase, mu 5 


NM_0 10360 


4E-51 




C3701 


Tnfa 


tumor necrosis factor-alpha 


NM_01 1659 


0.00009 


+ 


G0101 


Psme4 


proteasome activator subunit 4 


NM_134013 


4E-87 


+ 


G0701 


Mt2a 


metallothionein 2A 


NM_005953 


7E-20 




G1201 


Ggt1 


gamma-glutamyltransferase 1 


NM_013430 


3E-46 


+ 


G1501 


Txnrd2 


Thioredoxin reductase 2 


NM_006440 


9E-19 


+ 


G2001 


Pdia4 


Disulfite Isomerase 4 


NM_00491 1 


5E-46 


+ 


G2202 


Ahcy 


S-adenosylhomocystein hydrolase 


NM_0 16661 


1E-10 


- 


G2302 


Cyp2f2 


Cytochrome P450 2f2 


NM_007817 


9E-25 


- 


G3302 


Cytb 


Mesocricetus auratus cytochrome b 


YP_003208313 


9E-86 


++ 






4-Apoptosis 








A4301 


Pik3r1 


phosphatidylinositol 3-kinase 


NM_001 024955 


1E-80 


+ 


CI 701 


Tp53 


tumor supressor p53 


U07182 


6E-83 


- 


C7103 


Bel 10 


B-cell CLL/lymphoma 10 


NM_003921 


4E-37 


+ 


G3701 








5E-82 




G1503 


Nfkbl 


NF-kB 


NM_008689 


2E-102 


+ 


G2201 


Casp8 


Caspase 8 


NM_009812 


6E-29 


- 


G2305 


Eefld 


Eukaryotic translation elongation factor 1 delta 


NM_029663 


3E-17 


+ 


G3001 


Topors 


topoisomerase I p53-binding 


NM_1 34097 


4E-10 


- 


G4101 


Sh3kbp1 


SH3-domain kinase binding protein 1 


NM_031892 


1E-09 


+ 


G4102 








0.0008 




G5301 


Cmyc 


myelocytomatosis oncogene 


AJ582076 


2E-27 








5-Lipidogenesis 








A0501 


Pla2g2d 


phospholipase A2, group ID 


NM_01 1 109 


0.00003 


+ 


A3003 


Dhcr7 


dehydrocholesterol reductase 


NM_007856 


3E-36 


+ 


C0701 


Scl27a1 


solute carrier family 27A1 


NM_1 98580 


0.0 


+ 


C0901 


Acaal 


Acetyl-CoA acyltransferase 1 


NM_1 30864 


1E-77 


+ 


C4602 


Lpl 


lipoprotein lipase 


NM_008509 


1E-37 




G0201 


Star 


steroidogenic acute regulatory protein 


NM_0 11485 


7E-08 








6-Protein conformation or transport 








A3801 


Rpn1 


ribophorin I 


NM_1 33933 


3E-14 




A4701 


Ergic3 


ERGIC and golgi 3 


NM_1 98398 


6E-21 


+ 


A5302 


Gxyltl 


glycosyltransferase 8 domain containing 3 


NM_1 73601 


0.0002 


+ 


C0601 


Ppia 


peptidylprolyl isomerase A 


NM_021130 


7E-71 


+ 


C2001 


Grpell 


GrpE-like 1, mitochondrial 


NM_024478 


1E-13 


+ 
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Table 1 List of genes identified by Differential Display after DEHP exposure and classified by major biological func- 
tion, according to GO process (Continued) 



C2201 




1 1 dldAll 1 


MM OfWVM 


4E-33 




C3602 


Slc26a9 


col i itp r^rripr familv ?fS mpmhpr Q 

JUI U LC 1 1 ICI lull Illy 1 1 ICI 1 1 UCI J 


NM 177/43 

1 N J V 1 1 / / Z_^J 


1 E-9 




C7501 


Slc1 3a3 


cnliifp C7\xx\&x famik/ 13 mpmhpr 3 

jUIULC l_dl 1 ICI I CI 1 I 1 1 IV \3 IIICIIIUCI 3 


MM 

1 NIVI UJt"UJJ 


3E-29 




G1801 


Slc1 5a 1 


^nli itp r^rripr famik/ 1 ^ mpmhpr 1 
juiulc L-Ciiiici idiiiiiy i -J iiiciiiuci i 


MM 0S307Q 


0 0002 




G3702 


Rpn2 


rihnnhnrin II 

l l YJ\JyJ\ \\J\ 1 1 l ll 


NM 019642 


3E-26 




G6101 


Nrbpl 


nuclear receptor binding protein 1 


NM_013392 


4E-12 


+ 


G6601 


Slc6a8 


solute carrier family 6 member 8 

7-Cell cycle 


NM_005629 


0.0001 




A0401 


Ingap 


Mesocricetus auratus islet neogenesis associated protein 


U41738 


4E-79 


- 


A3002 


Cdkn2b 


cyclin-dependent kinase inhibitor 2B 


NM_004936 


1E-9 




A4801 


Ppplcc 


protein phosphatase 1, catalytic subunit, gamma 


NM_002710 


2E-78 


+ 


A5702 


Sec11a 


SEC1 1 homolog A 


NM_014300 


1E-19 


+ 


C4002 


Mapk4 


mitogen-activated protein kinase 4 


NM_002747 


9E-26 


+ 


C4701 


Ccndbpl 


cyclin D-type binding-protein 1 


NM_010761 


1E-17 


+ 


C6301 
C6302 


Zw10 


ZW10 homolog centromere/kinetochore protein 

O VJUICI lUllLUUIIb 


NM_004724 


1E-04 
1 E-05 






Rexo2 


nCAZ, nlNAA caUI lULIcdbc Z MUlllUIUy 


INIVI U 1 jjZj 


6E-21 




A5403 


Hk2 


1 IcAUM 1 Idoc Z 


mm nnniRQ 


4E-28 




A5602 


Pcbp2 


\jKj\y\\K-.) uiiiuiiiy piuLciii z 


MM 031QRQ 


1 E-34 




A7301 


Lyrm4 


1 YR mntif rnntaininn A 

l_ 1 1 X II IULII CUI ILall in iy " 


NM 201 3SR 

I NIVI Z.U I J JO 


0.0004 




C1801 


Itpripl2 


inositol 1,4,5-triphosphate receptor Interacting protein-like 2 


NM_001 033380 


4E-11 


+ 


C2401 


Tsn 


tr^nclin 

LI Ql \o\ 1 1 1 


NM 01 16S0 


4E-10 




1 


Trip4 


thi\/rf"iiH hir^rmr^no rarc±r\\T\r intorar-tr^r A 
LllylUlU lIUlIllUIlc IcLcpLUI IMLcldLLUI ^+ 


NM n \f0 13 

I nIVI U I UZ I 3 


n nnn^ 




C7102 


Tbrg3 


transforming growth factor beta regulated gene 3 


NR_027799 


0.001 


+ 


G0101 


Psme4 


proteasome activator subunit 4 


NM_014614 


4E-87 


+ 


G0501 


Gatad2a 


GATA zinc finger domain containing 2A 


NM_001 1 13346 


8E-15 


+ 


G0901 


HoxalO 


Homeobox A10 (HOXA10) 


NM_008263 


1E-47 


+ 


G1101 


lap 


Syrian hamster intracisternal A particle 


M10134 


3E-73 


+ 


G2303 


Zc3h12c 


zinc finger CCCH type containing 12C 


NM_001 162921 


7E-9 




G6801 


Fst 


follistatin 

Ribosomal proteins 


NM_006350 


3E-42 


- 


A3703 


Rpl37a 


ribosomal protein L37a 


NM_000998 


0.001 


+ 


A3901 
C3901 


Mrpl45 


mitochondrial ribosomal protein L45 


NM_025927 


2E-10 
6E-59 




A3902 


Rps21 


ribosomal protein S21 


NM_001024 


7E-14 


+ 


A4001 


./. 


mitochondrial 12S ribosomal RNA 


X84390 


0.0 


+ 


G2002 


Rpl27a 


60S ribosomal protein L27a 


NM_000990 


5E-18 


+ 


G3504 


RpHO 


ribosomal protein L10 


NR_026898 


4E-132 


+ 


G3601 
G3602 


Rpl28 


ribosomal protein L28 


NM_009081 


3E-22 
8E-180 


+ 


G5701 


Rpl22 


ribosomal protein L22 


NM_000983 


0.00009 


+ 



The DEHP effect is identified by (+) for 2.0-fold over-expression, (++) for 10-fold over-expression, (-) for 2.0-fold under-expression and (-) for 10.0-fold under- 
expression for at least one dose of DEHP. The 37 Differential Display fragments showing no match after comparison with the RefSeq database were not listed in 
the table. Id (Identification Number) represents the internal reference of the sequence used before the identification of Differential Display sequence. 



genes involved in promoter methylation {mettlS, 
mettSdl) were up-regulated. On the other hand, 7 genes 
were down-regulated (CCAAT/enhancer binding protein 
delta, methyl-CpG binding domain protein 1, Map 
kinase 3, Protease subunit 5, SWI/SNF related actin 



dependent regulator of chromatin, ATP dependent heli- 
case 4 and Forkhead A3). 

Xenobiotic metabolism genes such as cytochromes 
and glutathione S-transferases were also found to be dif- 
ferentially expressed, indicating a mobilization of cellular 
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defence and detoxication systems. An up-regulation of 
cyplbl and cyp2el was registered, whereas cyp2f2 was 
found to be down-regulated. Concerning GST, the Pi 
family was over-expressed while the Theta and Mu 
families were down-regulated. 

Differential Display results confirmed the down-regu- 
lation of c-myc and showed down-regulation of p53. A 
down- regulation of pro-apoptotic genes {casp8, topors..) 
and an over-expression of anti-apoptotic genes (bcllO, 
nfl<bl, sh3kbpl...) were also observed. 

Expression of genes involved in the regulation of the 
cytoskeleton by qPCR 

The mRNA level of the 21 genes involved in the regula- 
tion of the cytoskeleton that were identified as differen- 
tially expressed after 24 hrs was confirmed by qPCR. 
The expression of these genes was also studied after 
exposure to DEHP for 5 hrs. Out of the 21 genes, four 
(corolC, kif23, cdh3 and cttnbipl) were significantly up- 
regulated by DEHP treatment after 5 hrs of exposure 
(Figure 2) and one (nrp2) was significantly down-regu- 
lated. A clear dose-response relationship was observed 
for these 5 genes. After 24 hrs, these changes were con- 
firmed for 3 genes (up-regulation of corolC and kif23, 



and down-regulation of nrp2; Figure 3). Nevertheless, 
the down- and up-regulation was more pronounced 
after 24 hrs than after 5 hrs of DEHP exposure, for nrp2 
and kif23 respectively. For instance, in cells exposed to 
50 [iM of DEHP, Kif23 was up-regulated 17-fold at 24 
hrs versus 3-fold at 5 hrs. After 24 hrs, 5 other genes 
were significantly up-regulated (collAl, cripl, calm\3, 
dclk, and cttnbp2) by a factor ranging from 2.0 to 4.5 
with a dose-related effect (Figure 3). Eight other genes 
were significantly down-regulated (thbsl, flrt2, cdh3, 
has2, enah, ctnnbipl, lrrc8a and snx6), with an expres- 
sion ratio between 0.2 and 0.5 (corresponding to 5 and 
2 fold down-regulation respectively). All these genes 
were down-regulated in a dose-dependent manner, 
except for cdh3, enah, ctnnbipl, lrrc8a and snx6. A 
threshold was observed with the latter genes (12.5 [iM 
for cdh3; 50 [iM for enah, ctnnbipl, lrrc8a and snx6). 
Ctnnbipl was significantly down-regulated only for the 
lowest dose of DEHP (12.5 (iM). 

Although they had been identified as differentially 
expressed in DD, five genes (thyl, tubb2b, f5-actin, ple- 
khaS and nid2) were not shown to be significantly over- 
or under-expressed by qPCR. Yet the expression profiles 
of these genes indicated a dose-related increase for 




Control 

□ DEHP 0 |iM 

□ DEHP 12.5 \M 

■ DEHP 25 \iM 

■ DEHP 50 |iM 



Figure 2 Representative qPCR results of differentially-expressed genes involved in cytoskeleton regulation (according to the GO 
process), identified by Differential Display after 5 hrs of SHE cell exposure to DEHP. These histograms show the AACt score normalized by 
gapdh mRNA level. Error bars represent the standard deviation of the AACt score. Only mRNA levels showing a two-fold increase or decrease at 
least, were considered indicative (*) of a change in gene expression. 
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□ Control 

□ DEHPOuM 

□ DEHP 12.5 |iM 

■ DEHP25 nM 

■ DEHP50 nM 



Figure 3 Representative qPCR results of differentially-expressed genes involved in cytoskeleton regulation (according to the GO 
process), identified by Differential Display after 24 hrs of SHE cell exposure to DEHP. These histograms show the AACt score normalized 
by gapdh mRNA level. Error bars represent the standard deviation of the AACt score. Only mRNA levels showing a two-fold increase or decrease 
at least, were considered indicative (*) of a change in gene expression. 



tubb2b, fi-actin and pleckha5 but below the qPCR 2.0- 
fold threshold. As for thyl and nid2, the dose-related 
decrease was inferior to 0.5. 

Expression of apoptosis-related genes, PPARs and CYP4 
genes after DEHP treatment 

The expression level of bcl-2 and c-myc mRNA was 
used as controls of DEHP effects. An increased level of 
bcl-2 after 5 hrs of exposure and a decreased level of c- 
myc after 24 hrs (Figure 4) were observed according to 
qPCR, as expected. 

p53 was down-regulated in a dose- and time-depen- 
dent manner; a significant decrease of the mRNA level 
was found after 24 hrs at 50 [iM DEHP. 

None of the PPAR genes was identified as being differ- 
entially expressed by DD after DEHP exposure. In order 
to check these results, we measured the mRNA level of 
PPARa, PPAR p/S and PPAR 7, by qPCR using hamster 
specific primers. No change in the expression of these 
genes was observed by qPCR after 5 or 24 hrs of exposure 
with DEHP in our study conditions (Figure 4). The same 
verification was carried out for CYP4 genes. Neither Dif- 
ferential Display nor qPCR allowed us to identify signifi- 
cant expression changes compared to the control. 



Discussion 

The DDRT-PCR (Differential Display RT-PCR) techni- 
que was used in the present study to identify the differ- 
ential mRNA expression patterns between control and 
DEHP-treated SHE cells. Indeed, this technique is still a 
method of choice for non-sequenced or partially- 
sequenced organisms and is able to identify nonabun- 
dant, rare or novel transcripts [36]. 

Using Differential Display, we found 122 genes whose 
expression was altered by DEHP treatment (79 up-regu- 
lated and 43 down-regulated). The concentrations stu- 
died (12.5, 25 and 50 \iM) were in the range of 
concentrations that induced a morphological transfor- 
mation of SHE cells, i.e. concentrations up to 77 [iM for 
Mikalsen et al. [27] and in the range 25 (iM-150 [iM for 
Cruciani et al. [29]. 

We measured the mRNA level of genes involved in 
the regulation of the cytoskeleton using qPCR. This 
focus is justified by the fact that the modifications of 
cytoskeleton organization are early events in cell neo- 
plastic process [37] and can be recorded in SHE cells 
after 7 days of exposure to carcinogenic agents in cell 
transformation assays. Morphological transformation 
affects a few percentage of the mixed population of SHE 
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5 hrs 24 hrs 



Figure 4 Expression level of bcl-2, c-myc, p53, pi3kr1, PPARs and CYP4 family gene mRNA after treatment of SHE cells with DEHP, 
using qPCR. These histograms show the AACt score normalized by the gapdh mRNA level. Error bars represent the standard deviation of the 
AACt score. Only mRNA levels showing a two-fold increase or decrease at least, were considered indicative (*) of a change in the gene 
expression. We observed a significant increase in the bcl-2 mRNA level after 5 hrs of exposure and a significant decrease in c-myc and p53 mRNA 
levels after 24 hrs of DEHP exposition. PBkrl was found to be over-expressed for both lengths of exposure. None of PPARs or CYP4 genes was 
significantly over- or under-expressed after treatment. 



cells and all cell types [31]. From the present work, we 
can assume that the differentially-expressed genes mea- 
sured in the first 24 hrs of exposure reflect the first tar- 
gets of DEHP in the entire SHE cell population. The 
transcriptomic changes which were recorded correspond 
to the integrated mean of the cell responses significantly 
different in the exposed populations (p < 0.01), without 
consideration of cell specificity and sensitivity to DEHP. 
These significant expression changes in genes involved 
in cytoskeleton regulation, can be seen as early indica- 
tors of disturbances that will lead to cell transformation 
further in a few percentage of the most susceptible cells 
of the SHE population. The role of the cytoskeleton has 
been extensively studied in relation to invasion and 
metastasis, but little is known of its implication in the 
first stages of carcinogenesis. The identification of geno- 
mic changes associated with the triggering of cell trans- 
formation is useful from a mechanistic point of view 
and may be valuable in screening. 

Effects on cytoskeleton-related genes 

DEHP was shown to affect several functions related to 
the cytoskeleton. The genes involved in cytoskeleton 
regulation and identified by Differential Display are 
listed in table 2. To summarize, DEHP affects actin 



Table 2 List of genes involved in the regulation of the 
cytoskeleton and affected by DEHP 



5 hrs 



24 hrs 



Up-regulated genes corolC* (25 - 50 uM) 


corolC* (25 - 50 uM) 


kif23* (25 - 50 uM) 


kif23* (25 - 50 uM) 


cdh3* (12.5 - 25 - 50 


collar (50 uM) 


MM) 


cripr (25 - 50 uM) 


ctnnbipl* (25 uM) 


calml3* (12.5 - 25 - 50 




MM) 




dclkr (25 - 50 uM) 




cttnbp2* (50 mM) 




plekha5 




tubb2b 




p-octin 



Down-regulated 
genes 



nrp2* (25 - 50 piM) 



50 



nrp2* (25 - 50 uM) 
thbs* (50 uM) 
flrt2* (50 uM) 
cdh3* (12.5 - 25 - 
MM) 
hos2* (50 uM) 
enah* (50 mM) 
Ircc8a* (50 uM) 
snx6* (50 mM) 
ctnnbipl* (12.5 mM) 
nid2 
thyl 



This table summarizes the genes studied using qPCR. (*) indicates significant 
effects of DEHP (2.0-fold over- or under-expression) at concentration(s) 
specified in brackets. A trend for up- or down-regulation was found for the 
other genes (no concentration reported). 
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polymerisation and stabilization, as well as cell-to-cell 
and cell-to-matrix adhesion processes. The expression of 
genes involved in organelle transport, in cytoskeleton 
remodelling, or adhesion in response to external factors 
was also modified by DEHP. These results are in line 
with the recent findings of Posnack et al [38] who iden- 
tified disturbances in mechanical adhesion function and 
protein trafficking in rats cardiomyocytes exposed to 
DEHP. 

Actin polymerization and stabilization 

To summarize the basic process, actin polymerization 
requires the Arp2/3 complex that needs to be stabilized 
by Enable Homolog (Enah) and is regulated by coronins. 
Enah is involved in the dynamic reorganization of the 
actin cytoskeleton, and stimulates nucleation and poly- 
merization [39]. Coronins act on F-actin binding and 
bundling activities, but are able to inhibit the activity of 
Arp2/3 complex [40]. Actin polymers also require cor- 
tactin, which stabilizes nucleation sites for actin branch- 
ing and elongation [41,42]. Cripl facilitates actin 
filament bundling and stabilizes actin interaction with 
a-actinin too [43]. Linkage of actin polymers to 



adherens junctions, mainly composed of the transmem- 
brane proteins cadherins, is insured through binding to 
a-catenin and P-catenin [44]. 

Based on the gene expression data generated, we have 
tried to synthesize the effects of DEHP on actin organi- 
sation and cell adhesion specifically (Figure 5; over- 
expression in red; under-expression in green). A 5- and 
24-hrs exposure to DEHP over-expressed Coronin 1C 
(CorolC), resulting in F-actin disassembly [45]. Disorga- 
nization was amplified by under-expression of Enah 
involved in actin nucleation and polymerization, and 
expression of Cttnbp2 that counteracts cortactin which 
is known to stabilize the actin network. On the other 
hand, the binding of actin filaments to cadherins 
through catenin links appears to be reinforced owing to 
under-expression of Ctnnbipl (a P-catenin blocker) and 
over-expression of Cripl, which intensifies fixation to 
actinin. Globally, the effects of DEHP on actin cytoske- 
leton disturb actin polymerization while intensifying 
binding on actinin and catenins. Posnack et al. [38] 
explored DEHP effects on rats cardiomyocytes in a 
range of concentrations two and three orders of 



Response to external factors 




protein which gene have 
not be identified in 

Differential Display Extra Cellular Matrix 

■ overexpression 

■ underexpression 

Figure 5 Representative scheme of the genes affected by DEHP and involved in cytoskeleton regulation. Genes in red have been found 
to be over-expressed using Differential Display (DD) and genes in green under-expressed. Dotted lines represent genes identified by DD but 
whose expression was not significant in qPCR. 
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magnitude higher than here. They found an over- 
expression of actinin, a-catenin and N-cadherin in a 
concentration-dependent manner. 
Cell-cell and cell-matrix adhesion 

Cell-cell adhesion and cell-matrix adhesion were also 
affected by DEHP treatment. The decrease in the P- 
Cadherin (Cdh3) mRNA level after 24 hrs of exposure 
indicates that DEHP weakened cell-cell contact, after a 
transient increase at 5 hrs of exposure for all doses 
tested. Weakening of cell-matrix adhesion may result 
from a decrease in the Hyaluronan synthase 2 (has2) 
mRNA level [46] and in Thrombospondin (Thbsl), an 
adhesive protein that interacts with fibronectin, laminin, 
integrins and collagen [47]. Loss of cell adhesion may 
also be explained by over-expression of CorolC because 
this gene negatively regulates cell-matrix adhesion 
through focal adhesion kinase (FAK) -mediated signalling 
[45]. Also, under-expression of Enah, which is known to 
be involved in the control of cellular adhesion by the 
recruitment of proteins containing SH3- domain [48], 
contributes to the loss of cell-cell adhesion. 

In addition, DEHP may lessen extracellular matrix 
adhesion by reducing the expression level of a number 
of transmembrane proteins involved in cell-matrix con- 
nections: Fibronectin leucine rich 2 (Flrt2) and Leucine 
rich repeat 8A (Lrrc8a) [49], Nidogen 2, which connects 
laminin- 1 to the matrix [50], and Thy-1, which mediates 
fibroblastic adhesion [51] and is Thbsl expression- 
dependent [52]. On the other hand, DEHP effects rein- 
force the extra-cellular matrix through an over-expres- 
sion of collAl increasing collagen. This effect may be 
seen as a compensatory reaction to the weakening of 
cell-to-matrix link proteins by DEHP. Sobarzo et al. [53] 
demonstrated an up-regulation of N-cadherin and a- 
catenin in rat testis after 2 and 7 days of DEHP treat- 
ment, suggesting also a deregulation of cell adhesion 
molecules in seminiferous tubules. 

DEHP decreases the response to external factors, 
such as the Vascular Endothelial Growth Factor 
(VEGF) or the Epidermal Growth Factor (EGF) 
through under-expression of neuropilin 2 (nrp2) and 
sorting nexin 6 (snx6) respectively. Nrp2 is a mem- 
brane receptor capable of binding VEGF and sema- 
phorins, therefore its under-expression may inhibit cell 
adhesion and migration via the loss of integrins [54]. 
Snx6 is able to interact with EGF receptor and Trans- 
forming Growth Factor (TGF)-P receptor [55]. Under- 
expression of snx6 and thbsl may lead to decreased 
interaction with Latent TGF Binding Protein (LTBP) in 
the upstream of the TGF-P pathway contributing to 
the repression of the TGF-P signaling pathway [56,57]. 
Under-expression of TGF-P is known to decrease 
apoptosis in rodent hepatocytes treated with peroxi- 
some proliferators (PPs) [24]. 



Organelle transport and cytoskeleton remodelling 

DEHP also interferes with functions of microtubules 
(composed of a- and p -tubulin). Kif23, which encodes a 
kinesin protein, was highly over-expressed after 5 hrs 
and 24 hrs of DEHP exposure. Kif23 has been shown to 
transport membranous organelles and protein com- 
plexes from cell nucleus to cell periphery in a microtu- 
bule- and ATP-dependent manner. Doublecortin-like 
kinase (Dclk) is a microtubule-associated protein encod- 
ing a Ca 2+ /calmodulin-dependent kinase. Its activities 
on binding and microtubule polymerization facilitate 
cell motility by remodelling the microtubule cytoskele- 
ton [58]. Over-expression of dclk at 24 hrs of DEHP 
treatment is in line with an increased trend in p -tubulin 
{tubb2b). 

Calmoduline-like 3 (calml3) was over-expressed after 
24 hrs of DEHP exposure. Calmodulin (CaM) is a cal- 
cium-binding protein that translates the Ca 2+ signal into 
a wide variety of cellular processes, including the regula- 
tion of cytoskeleton remodelling acting with Caldesmon 
[59] or with Wnt pathway [60]. Calml3 is a CaM family 
member protein which increases cell motility by stabiliz- 
ing and increasing myosin- 10 for cell migration [61]. 
Other genes involved in signal transduction pathways and 
cytoskeleton regulation 

We measured an over-expression level of phosphatidyli- 
nositol 3-kinase rl (pi3krl) using Differential Display 
and qPCR. Pi3k is a key signalling molecule in the PIP 3 
signalling transduction pathway and in actin reorganiza- 
tion and cell adhesion [62] and is able to regulate the 
synthesis of collagen I [63]. An activation of PI3K is also 
associated with a phosphorylation-dependent activation 
of Akt which contributes to tumorigenesis and metasta- 
sis [39]. The over-expression of pi3krl can be related to 
the under-expression of ctnnbipl which interacts with 
P-catenin. In addition to the function of P-catenin in 
the actin cytoskeleton, its role in the regulation of Akt 
pathway activation [64] or in Wnt pathway regulation 
[65] is advanced. This protein forms part of a complex 
that captures growth and proliferation signals from the 
cell surface and is then activated to stimulate the 
expression of genes involved in cell proliferation. It 
would be worth studying P-catenin-dependent transcrip- 
tion in relation to carcinogenicity. 

DEHP effects in the SHE model compared to rats and 
mice 

While the expression of cyplbl and cyp2el was up- 
regulated and cyp2f2 under-expressed, no change in 
expression level of CYP4 genes was found using DD and 
qPCR after DEHP exposure in our experimental condi- 
tions. CYP4 genes are said to be involved in peroxisome 
proliferation [66]. Eveillard et al. [67] who studied the 
involvement of DEHP in lipidogenesis in rats, found a 
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slight increase in the PPARa level after 21 days of oral 
exposure to DEHP (200 mg/kg.day). They registered a 
significant increase in CYP4 levels after 14 days [68] and 
after 21 days [67] of exposure. On the other hand, we 
found no increased mRNA level of CYP4 and PPAR 
genes in DEHP-treated SHE cells. This underlines that 
the genes expression changes noted in the present study 
are independent of PPARs induction. Eveillard et al. [67] 
found that induced expression of cyp2bl0 by DEHP was 
also independent of PPARa induction but CAR depen- 
dent. No change in CAR expression was registered in 
SHE cells, which may explain why no change in cyp2bl0 
was noted. Our results are consistent with the study of 
Ren et al. [20] who identified DEHP regulated genes 
independent of PPARa and CAR in rats and mice. 

In our study, lipogenesis and xenobiotic metabolism 
pathways were impacted by DEHP, but not in a major 
prior way. This may be explained by the lower sensitivity 
of the hamster model compared with rats and mice to 
peroxisome proliferators [69]. Indeed, the Syrian hamster 
model presents an intermediate response between rats or 
mice and humans who are known to be non-responsive 
to PP induction [24]. The hamster model, like humans, is 
less responsive to PP induction than rats and mice, which 
is an advantage for mechanistic studies of PP effects and 
for screening human chemical carcinogens. 

On the other hand, three genes (Actin p, Lipoprotein 
lipase, and Acetyl-Coenzyme A acyltransferase la) and 5 
gene isoforms (Glutathione S-transferase-pi, -mu, -theta, 
CCAAT/enhancer binding protein and Nuclear receptor 
subfamily 1) were commonly found in our study and 
those carried out by Eveillard et al. [67,68], suggesting a 
pattern of response specific to DEHP. 

Takashima et al. [70] also found similar responses in 
DEHP-treated mice. Up-regulation of rablb, a RAS 
oncogene family member involved in cellular signal 
transduction or survival, was found in the latter study 
and in the present one. P -Tubulin was clearly over- 
expressed in mice, a trend which was noted in our 
study. Some gene isoforms of cadherin, nidogen, cypl 
family genes or LIM domain were also impacted in the 
liver of mice exposed to DEHP [70]. 

DEHP effects on transcription factors 

Other genes identified by Differential Display and 
involved in transcription and signal transduction path- 
ways or apoptosis were also targeted by DEHP. A signif- 
icant under-expression of pS3 was found after 24 hrs of 
DEHP exposure using Differential Display and qPCR. 
This under-expression is in line with the anti- apoptotic 
effects of DEHP. 

We confirmed the over-expression of bcl-2 after 5 hrs 
and the under-expression of c-myc after 24 hrs, events 
reported in a previous study on DEHP treated SHE cells 



in conditions similar to the present ones [26]. Map 
Kinases such as Mapk3, Mapk4 and Mapkl5 were tar- 
geted by DEHP. Further investigations of Map Kinase 
pathways could be relevant due to their involvement in 
activities of transcription factors. 

The G protein-coupled estrogen receptor (gper) was 
found to be over-expressed in Differential Display. Gper 
can be activated by estrogen-like compounds and its 
effect on cytoskeleton architecture has been reported 
[71]. Because of its implication in the regulation of 
MAPK [72] or TGF-p pathways [73], it would be worth- 
while to investigate gper further. 

Performances of DD 

The confirmation of differentially-expressed genes by 
qPCR showed that the expression levels of more than 
75% of genes identified by DD were confirmed by qPCR. 
A comparative table of the sensitivity of DD versus qPCR 
is given in additional file 1. qPCR is more likely to quan- 
tify subtle changes in the expression level of mRNAs at 
different concentrations while DD seems to be more sen- 
sitive but is less discriminating. To summarize, 35% of 
the genes identified as differentially expressed in DD gave 
the same response at the same DEHP concentrations 
with qPCR while 40% were detected by DD at a lower 
DEHP concentration than with qPCR. 

Conclusion 

Transcriptional responses of SHE cells to DEHP were stu- 
died in conditions inducing the cell neoplastic transforma- 
tion, in order to identify gene expression changes in 
relation with effects of this non-genotoxic carcinogen. 
Functions impacted by DEHP were found to be PPAR- 
independent. Effects on cytoskeleton related genes indi- 
cated disturbances on actin polymerization and stabiliza- 
tion, cell-cell and cell-matrix adhesion and protein 
trafficking. 

This is the first study that elucidates the genomic 
changes of DEHP on the organization of the cytoskele- 
ton. Whether the expression changes of cytoskeleton- 
related genes identified here such as corolC, nrp2, kif23, 
are specific to DEHP or to cell transforming agents more 
generally would require further studies. To answer, the 
gene sets identified as significantly over- or under- 
expressed in this study must be explored on other non- 
genotoxic carcinogens to identify biomarkers predictive 
of early events in the multistep carcinogenic process. 
Early disturbances in the expression of cytoskeleton- 
related genes should be considered good candidates. 

Methods 

Chemicals 

DEHP (C.A.S. No. 117-81-7, purity 99%), purchased 
from Aldrich Chemicals (Gillingham, England) was 
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dissolved in the DMSO solvent (C.A.S. No. 67.68.50). 
The latter was obtained from Sigma Aldrich (St Quentin 
Fallavier, France) and was used at a final concentration 
of 0.1%. 

Nucleic acid stain Gelred purchased from Interchim 
(Montlucon, France) was used at a final concentration 
of 1:10000 (v/v). 

All chemicals used for this study were electrophoresis 
grade or molecular biology grade. Their origin is speci- 
fied in the following sections. 

SHE cell culture and treatment 

SHE cells were isolated from Syrian hamster embryos at 
day 13 of gestation using the procedure described by 
Pienta et al. [74] and in accordance with the modifica- 
tions suggested by Elias et al. [75]. Differentiated tissues, 
such as eyes, heart and viscera were removed and the 
remaining tissues were dissociated by dispase (1.2 U/ 
ml). Stock cells were preserved in liquid nitrogen. 

After thawing, passage-2 primary cells were pre-cul- 
tured until they reached 80% confluency. The culture 
medium was Dulbecco's modified Eagle medium 
(DMEM) (Gibco, Invitrogen, Cergy Pontoise, France) 
supplemented with 10% foetal calf serum (FCS) 
(Hyclone, Brebieres, France; lot N°#ASB28835), 1.5 g/L 
NaHC0 3 at 37°C in a 10% C0 2 humidified atmosphere 
and pH 7.0. No phenol red was added to the medium. 

Cells used for the DEHP studies were sampled from a 
monolayer during the growing phase, 48 hrs after seed- 
ing. Cells were trypsinized and treated during replating 
with DEHP at concentrations of 0 (iM (vehicle control), 
12.5 (iM, 25 [iM and 50 |iM in DMEM culture medium 
supplemented with 10% FCS. Cells were then incubated 
for 5 hrs and 24 hrs at 37°C in a 10% C0 2 humidified 
atmosphere. 

RNA isolation 

Total RNA extractions were performed directly in the 
dish, using Nucleospin RNA II Extract Kit (Macherey 
Nagel, Hoerdt, France), according to the manufacturer's 
instructions. A DNAse I treatment was performed 
directly through the column used to collect RNAs and 
before the elution phase of DNA-free RNA. 

RNA was quantified by spectrophotometry (Nanodrop, 
Labtech) measuring the A 260 /A 280 ratio and its quality 
was ensured by electrophoresis using a 1% RNase-free 
agarose gel. Aliquots were stored at -80°C before use for 
Differential Display and Real-time PCR. 

Anchored Reverse transcription (RT) and Differential 
Display 

The Differential Display was performed as described by 
Liang et al. [35], with minor modifications concerning 
DD fragment revelation with GelRed. 



For Differential Display, three separate RT reactions 
were performed with a different one-base anchored 
oligo-dT primer (H-dTA, H-dTC and H-dTG) to pro- 
duce three different subsets of cDNA pools. The 
sequences of the anchored and the arbitrary primers are 
given as additional file 2. The RT reactions were carried 
out using 2 uL of each primer (50 |iM) and 4 |ig of total 
RNA. 8 (iL of RevertAid M-MuLV Reverse Transcrip- 
tase 5x reaction buffer (Fermentas, Saint-Remy-les- 
Chevreuse, France), 1.5 (iL of 10 mM dNTPs (Fermen- 
tas, Saint-Remy-les-Chevreuse, France) and up to 35 |iL 
Nuclease-free water were added to each tube, mixed, 
then heated at 70°C for 3 min. Tubes were centrifuged 
and incubated on ice for 5 min, then 2 \iL (40 U) of 
RNaseOUT Recombinant RNase Inhibitor (Invitrogen, 
Cergy Pontoise, France), 1 (iL (200 U) of RevertAid M- 
MuLV RT (Fermentas, Saint-Remy-les-Chevreuse, 
France) and 2 (iL of Nuclease-free water were added to 
each tube. Each tube was mixed by gentle pipetting then 
incubated in a thermocycler at 42°C for 1 h, followed by 
95°C for 10 min. The tubes were then centrifuged and 
stored at -80°C until use. 

Amplification was then performed using combinations 
of the three original anchored primers from the reverse 
transcription step and eighty arbitrary 13-mers (H-AP), 
giving a total of 240 amplification combinations. All 
reactions contained 2 |iL of a lOx PCR buffer containing 
25 mM of MgCl 2 , 1.6 |iL of 1 mM dNTP mix, 1 U of 
Taq Polymerase (Euromedex, Souffelweyersheim, 
France) and the primer combination at a final concen- 
tration of 1 (iM. Tubes were incubated for 5 min at 95° 
C. The next 40 cycles were 95°C for 30 s (denaturation), 
40°C for 2 min (annealing), 72°C for 1 min (amplifica- 
tion). A final extension of 72°C for 10 min completed 
the cycle. 

After thermocycling, PCR-amplified fragments were 
resolved in a 6% native polyacrylamide gel in 1 x TBE 
buffer (89 mM Tris base/89 mM boric acid/2 mM 
EDTA, pH 8.0), using 10 |iL of PCR product mixed with 
2 |il of loading buffer (0.05% xylene cyanol, 40% sucrose, 
20 mM EDTA, pH 8.0). Gels were run at 100 V for 20 
hrs, then the fragments were stained with GelRed IX in 
water, for 30 min in the dark. Bands on the gel were 
revealed on a UV-transluminator. PCR products that 
showed differential expression between control and trea- 
ted samples were identified with QuantityOne® 1-D 
analysis software (BioRad, Marne-la-Coquette, France). 

Bands which were up- or down- regulated more than 
2-fold, were selected and characterized in the next step 
of analysis. Differentially-expressed bands were excised, 
reamplified and their sizes were checked before cloning. 
To summarize, fragments of interest were recovered 
using a clean razor blade and extracted from the gel 
matrix by boiling in 200 (iL of buffer (10 mM Tris/1 
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Table 3 List of the primers used for real-time qPCR 



Genes 



Accession N° 



Primers sequence (5'-3') 



Housekeeping gene 
gapdh 



Glyceraldehyde 3-phosphate dehydrogenase 



Regulation of cytoskeleton 
tubb2b 

thyl 

actin 

collal 

thbsl 

plekha5 

kif23 

has2 

flrt2 

enah 

dclkl 

ctnnbipl 

cripl 

co role 

cdh3 

calml3 

cttnbp2 

Irrc8a 

snx6 

nrp2 

nid2 

Regulation of apoptosis 
pi3kr1 

bcl-2 

c-myc 

p53 



p-Tubulin 
Thy1 Antigen 

p-Actin 
Collagen a1 
Thrombospondin 1 
Plekstrin homology A5 

Kinesin 23 
Hyaluronan synthase 2 
Fibronectin leucine rich 2 
Enabled homolog 
Doublecortin like 1 
Catenin |3 interacting protein 1 
Cystein rich protein 
Coronin 
Cadherin 3 
Calmodulin like 3 
Cortactin binding protein 2 
Leucine rich repeat 8A 
Sorting nexin 6 
Neuropilin 2 
Nidogen 2 

phosphatidylinositol 3-kinase (p85a) 

B-cell lymphoma 2 
cellular myelocytomatosis oncogene 
tumor supressor p53 



PPARs genes 
ppar a 



peroxisome proliferator activated receptor alpha 



DQ403055 

NM_023716 

NM_009382.3 
NM_007393 
NM_000088 
NM_011580 
NM_019012 

NM_024245.4 
NM_008216 
NM_201518 
NM_008680 
NM_019978 
NM_020248 
NM_007763 
NM_0 14325 
NM_007665 

NM_027416.3 
NM_030249 
NM_1 77725 
NM_021249 

NM_00 10774 
NM_008695 

NM_001 077495 
AJ582074 
AJ582076 
U07182 

AY 170844 



F: CAATGACCCCTCATOACC 
R: GACAAGOTCCCGTCTCAG 

F: GCAACATGAATGACCTGGTG 
R: ACCAGAGACCCAGCACAAAC 

F: AAGGCCTCTGCCTGTAGTGA 
R: GAAGAGGCAGGTOCAAGAC 

F: CACCACCACAGCCGAGAG 
R: CCAGGGAGGAAGAGGATGC 

F: GGGTCAmCCACATGCm 
R: TCCGGGmCAGAGTACCAC 

F: CCAAAGCCTGCAAGAAAGAC 
R: CCTGOTGTCCAAAOTGA 

F: GTGCATCTGCCTGAAGACAA 
R: TGGGAACCmAACGACTGG 

F: CCTGAGCmCCTGACCAAG 
R: AGTCOTCTGGGTGGTGTG 

F: CGGAGGACGAGTCTATGAGC 
R: ^CCGGTGTOCAAAAAG 

F: ACCGCACTGTGGAAGATACC 
R: GCAAGACAACGAGCACAAAA 

F: GCCTATGOTCAGCAOTCC 
R: GGGCGATOTOTCTGACAT 

F: AGCCTCCACCAGCTCAG™ 
R: CCATACACATCGCTCCATO 

F: TOGCTGCAGAAAGAAACCT 
R: CCAGCCAATCACAACC^ 

F: AGTCCAGAGCCTGCAACCTA 
R: GGAGTAGCAGGGATGATOC 

F: GCAGAAGAGTGGTOGAAGG 
R: TGATCAGGTCGCAOTOTG 

F: CACACGACCTCATGTOACC 
R: CTGTACCTCATGGCCCAOT 

F: ATCGACAAGGATGGAAACG 
R: ATCTACCTCCTCGTCGCTCA 

F: GACAAAGAAGGCTGGACTGC 
R: CTCACCCACGGAAATCC™ 

F: AGAGCCCACTOCCCCAACT 
R: GTGTGCAGAAGCACGAGGTA 

F: CCAAGACCTGA^GATGOTC 
R: CATGCATCGCAACTGTCTO 

F: ATAAGCACTGATGTCCCACTG 
R: GAGTOCTCCAATCTCOTCA 

F: CTCATOAGTOTGCCTGC 
R: ATAGCTGCCTCATGACATCG 

F: CGAGCCCGACCGGAGGTGAA 
R: CGCACACTGCCGTCCGAG^ 

F: CGCAGAGATGTCCAGTCAGC 
R: CGAACTCAAAGAAGACCACAA 

F: GACCCTGATCGGACCTOT 
R: CGACTCCACAGCOTCTCTC 

F: ATGACGGAAGTOTAAGAC 
R: TCGGATAAGATGCTGAGG 



F: GmcmCGGCGAACTA^ 
R: ACACGTGAGAATCTCTGGT 
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Table 3 List of the primers used for real-time qPCR (Continued) 


ppar p/S peroxisome proliferator activated receptor beta/delta 


AF486582 


F: 


TGCAAGATCCAGAAGAAGAA 






R: 


GTAGATGTGGTGGAGAAGG 


ppar y peroxisome proliferator activated receptor gamma 


AB525757 


F: 


GGACCTCTCTATGATGGATG 






R: 


GGATGCAGGTOTACmGA 


CYP4 family genes 








cyp4a17 Cytochrome P450 4A17 


AJ555628 


F: 


ACCAGATGCCCTACACTACC 






R: 


GTGCGTAAATGGAGAGTACA 


cyp4a!8 Cytochrome P450 4A18 


AJ555629 


F: 


ACCAGATGCCCTACACTACC 






R: 


GGCCATAAATGGAGATOCA 


cyp4a19 Cytochrome P450 4A19 


AJ555630 


F: 


ACCAGATGCCCTACACTACC 






R: 


GTGCATAAATGGAGAGTGTG 



mM EDTA/1% SDS (w/v), pH 8.0) for 15 min. After 
overnight precipitation at -80°C, the eluted DNA was 
reamplified using the same primers and PCR conditions 
as the ones used in the DD-PCR step. Reamplified DNA 
was run in a 1.5% agarose gel containing IX GelRed 
and recovered using NucleoSpin ® Extract II kit 
(Macherey Nagel, Hoerdt, France) before cloning. 

Cloning was carried out using a TA Cloning Kit 
(pGEM-T, Promega, Charbonnieres, France), according 
to the manufacturers instructions. Plasmid DNA was 
extracted from the cultures using Nucleospin ® Plasmid 
QuickPure (Macherey Nagel, Hoerdt, France), according 
to the manufacturers instructions and sequenced bidir- 
ectionally by the DNA sequencing service of MWG 
Operon (Ebersberg, Germany), using T7 and SP6 
primers. 

Identification of differentially-expressed genes 

Sequences were compared with the National Centre of 
Biotechnology Information Gene Bank database (http:// 
www.ncbi.nlm.nih.gov) using the tBLASTx algorithm 
and RefSeq mouse or Refseq human as a reference. 

Confirmation of differentially-expressed sequences by 
Quantitative Real Time PCR (qPCR) 

First-strand cDNA was synthesized from 2 ug of total 
RNA using VILO-Superscript™ III reverse transcriptase 
(Invitrogen, Cergy Pontoise, France) and random-hex- 
amer primers. To summarize, 2 ug of total RNA was 
combined with 4 uL of 5X VILO™ reaction mix (con- 
taining RT buffer, MgCl 2 , dNTPs and random primers) 
and 2 uL of 10X enzyme mix (containing Superscript® 
III and RNase inhinitor). The final volume was adjusted 
to 20 uL and the reaction mix was incubated at 42°C for 
60 min. Then, cDNAs were diluted 20-fold, according to 
the manufacturer's instructions, before qPCR amplifica- 
tions. The oligonucleotides used as primers in the quan- 
titative real time PCR assay are described in table 3. If 
possible, at least one primer in each pair spanned an 
exon-intron boundary. PCR was carried out using Fast 
SYBR®Green Master Mix (Applied Biosystem, 



Courtaboeuf, France). Amplifications were performed on 
a StepOnePlus Real-Time PCR system (Applied Biosys- 
tem, Courtaboeuf, France). Each qPCR reaction con- 
tained 10 uL of 2X Fast SYBR®Green Master Mix, 5 uL 
of primers, 2 uL of diluted cDNA and 3 uL of Nuclease- 
free water. Amplification parameters were set as follows: 
initial denaturation (95°C, 3 min), and then amplifica- 
tion (95°C, 3 s and 60°C 30s) for 40 cycles. Glyceralde- 
hyde 3-phosphate dehydrogenase (gapdh) mRNA level 
was used as a housekeeping gene to normalize qPCR 
data. This gene was chosen because DEHP exposure did 
not affect its expression unlike P-actin which was also 
tested (data not shown). qPCR results were analyzed 
using the software provided with the thermocycler and 
DataAssist, using the AACt method [76]. Each validated 
primer pair used yielded a single peak of dissociation on 
the melting curve. The efficiency calculated by standard 
curve with five log- 10 dilution points was between 0.95 
and 1.05. A 2.0-fold threshold and a p-value of 0.05 
were used to determine the significance of differential 
expression levels according to the standard parameters 
of DataAssist. 

Additional material 



Additional file 1: Comparative table of the sensitivity of DD versus 
qPCR. Comparisons between the ratio of bands intensity compared to 
the control on the gel using an image analysis software (QuantityOne®® 
1-D analysis software BioRad, Marne-la-Coquette, France) for Differential 
Display and the AACt score normalized by gapdh mRNA level after 
analysis with StepOne and DataAssist (Roche Applied Biosystem, 
Courtaboeuf, France) for qPCR. 

Additional file 2: Primers for Differential Display. List of the 
sequences of the anchored and the arbitrary primers used for the 
Differential Display experiments. 
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